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Abstract 

Dengue fever is among the most widespread vector-borne infectious diseases. The primary vector of dengue is the Aedes 
aegypti mosquito. Ae. aegypti is prevalent in the tropics and sub-tropics and is closely associated with human habitats 
outside its native range of Africa. While long established in the southeastern United States of America where dengue is re- 
emerging, breeding populations have never been reported from California until the summer of 2013. Using 12 highly 
variable microsatellite loci and a database of reference populations, we have determined that the likely source of the 
California introduction is the southeastern United States, ruling out introductions from abroad, from the geographically 
closer Arizona or northern Mexico populations, or an accidental release from a research laboratory. The power to identify 
the origin of new introductions of invasive vectors of human disease relies heavily on the availability of a panel of reference 
populations. Our work demonstrates the importance of generating extensive reference databases of genetically 
fingerprinted human-disease vector populations to aid public health efforts to prevent the introduction and spread of 
vector-borne diseases. 
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Introduction 

Dengue fever is re-emerging in the United States of America 
(USA) after many years of absence [1], a trend also observed 
around the world [2] . Risk for dengue infection coincides with the 
distribution of mosquitoes capable of transmitting dengue virus 
(DENV). In most areas of the world, Aedes aegypti is the mosquito 
species responsible for DENV transmission. Ae. aegypti is a 
domestic species in the sense that, outside Africa, it is closely 
associated with human habitats and is often transported by 
humans and their commerce [3]. In the USA, Ae. aegypti is 
currently found in the southern states, Arizona eastward, with 
year-round breeding confined to latitudes below 33°N [4]. The 
state of California has an active and extensive mosquito- 
monitoring program since 1917 and in the past has only detected 
sporadic specimens of Ae. aegypti near airports [5]. Confirmed 
breeding populations of Ae. aegypti in California were never 
reported until the summer of 2013, when they were found in the 
Central Valley counties of Fresno and Madera and the coastal 
county of San Mateo (Figure 1C). When exotic vectors of 
important diseases are introduced into a new and densely 
populated region like California, it is important to understand 
where they originated in order to implement appropriate 
containment, control, and elimination strategies. As part of our 
ongoing studies of the genetic variation o£Ae. aegypti populations, 
we report here the use of 12 highly polymorphic microsatellites 
and a database of reference populations to identify the likely origin 
of the California Ae. aegypti detected in 2013. The power to 



identify the source of the introduction relies heavily on the 
completeness of the reference panel and thus, the origin could be 
further narrowed down as the reference database gets expanded. 

Methods 

Aedes aegypti collections 

Aedes aegypti adults and larvae for this study came from 
collections made between 2004 and 2013 from 33 locations across 
1 1 countries in the Americas and Asia (Table SI). Mosquitoes that 
were genotyped arrived as either eggs from oviposition traps, or as 
larvae or adults in 70-100% ethanol. Eggs were hatched at the 
Yale School of Epidemiology and Public Health insectary and 
reared to adults for identification and preservation in 100% 
ethanol at — 20°C. Mosquitoes included in this study came directly 
from the field, from multiple oviposition traps, unless indicated in 
Table SI. 

DNA extraction and genotyping 

Total nucleic acids were extracted from 1205 individual Aedes 
aegypti mosquitoes with the DNeasy Blood and Tissue kit 
(Qiagen), according to manufacturer instructions. Samples were 
treated with 4 \ll of RNAse A (Qiagen) and stored at -20°C until 
further analysis. Our reference dataset included eleven popula- 
tions of Ae. aegypti from around the world that have been 
screened for allele frequencies at 12 highly polymorphic 
microsatellite loci in previous studies [6,7]. Additionally, we 
genotyped 1 9 natural populations from Asia and the Americas as 
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Author Summary 

Dengue fever is a vector-borne disease that has dramat- 
ically increased in frequency and severity in recent 
decades, affecting over half of the world's population. 
The primary vector of dengue is Aedes aegypti. Breeding 
populations of this mosquito were reported for the first 
time in California in summer 2013, and the origin of these 
populations was unknown. We genotyped Ae. aegypti from 
California and compared their genetic composition and 
diversity against reference populations from Asia and the 
Americas, and with laboratory strains, to identify the 
source of this introduction. We identified the southeastern 
United States of America as the likeliest origin of California 
Ae. aegypti. This information will aid the design of 
containment, control, and elimination strategies to pre- 
vent the establishment of a permanent breeding vector 
population in the state. Note that the resolution achieved 
in tracking the origin of new disease-vector introductions 
is proportional to the completeness of the reference 
dataset. Our analysis was delayed to gather further data on 
regions close to California and would otherwise have been 
completed within a week. This work highlights the 
importance of having an extended dataset of human 
disease-vector genotypes in place prior to the occurrence 
of invasions that compromise public health. 

well as samples of 3 populations from California: Madera, Fresno, 
and San Mateo; for the same 12 loci (Table SI). Individual 
mosquitoes were genotyped as described in Brown et al. [6] . The 
microsatellite loci analyzed were: Al, AC2, CT2, AG2, B2, B3, 
A9, AC 4, AC1, ACS, AG1, and AG4 (Table S2). Briefly, 
polymerase chain reactions were conducted as 10 ul reactions 
using the Type-it Microsatellite PCR Master Mix (Qiagen), 
25 nM of each forward primer, 250 nM of each reverse primer, 
and 500 nM of a fluorescently labeled Ml 3 primer. Thermo- 
cycler conditions were as follows: 94°CxlO', 35 x(94°C x30", 
54°Cx30", 72°Cx30"), and 72°Cx5'. The resulting products 
were processed for fragment analysis at the DNA Analysis Facility 
at Science Hill at Yale University, using GS 500 Rox internal size 
standard (Applied Biosystems). Microsatellite alleles were scored 
using GeneMapper v4.0 (Applied Biosystems). 

Analyses 

All microsatellite loci were analyzed for within-population 
deviations from Hardy-Weinberg equilibrium (HWE) by F is , as 
well as for linkage disequilibrium (LD) among loci pairs, using the 
online version of the GENEPOP software [8,9] with 10,000 
dememorizations, 1,000 batches, and 10,000 iterations per batch 
for both tests. The G; s , analogue to F; s , was also calculated in 
GenoDive v.2.0b25 [10] to evaluate deviations from HWE using 
10,000 permutations. A Bonferroni correction was applied to the 
resulting matrix to correct for multiple testing. The Adegenet 
package v. 1.3.9. [11] available for the R software v. 3.0.1. [12] 
was used to compute the number of alleles at each microsatellite 
loci and the average observed and expected heterozygosities per 
population. Principal Component Analysis was implemented with 
the same package. Allelic richness and private allelic richness were 
calculated in HPRARE [13], which uses rarefaction to correct for 
unequal sample sizes. 

Population structure and assignments of individuals from the 
three California populations to their genetic cluster of origin was 
performed via the Bayesian clustering method implemented by the 
software STRUCTURE v. 2.3 [14]. STRUCTURE identifies 
genetic clusters and assigns individuals to these clusters with no a 



priori information of sample location. The most likely number of 
clusters (K) was determined by conducting 20 independent runs 
from K = 1 to 29 on individuals of all American and Asian 
populations, and from K = 1 to 2 1 on individuals of the North 
American cluster. Each run assumed an admixture model and 
correlated allele frequencies using a burn-in value of 100,000 
iterations followed by 500,000 repetitions. The optimal number of 
K clusters was determined following the guidelines of Prichard 
et al. [14] and the Delta K method from Evanno et al. [15,16]. 
Results were plotted with the program DISTRUCT v. 1.1 [17]. 

GENECLASS2 [18] was used to perform individual and group 
assignment tests on the three California populations against the 
reference population dataset using the Bayesian criteria for 
likelihood estimation [19]. The reference dataset included all 
populations identified in STRUCTURE as belonging to the same 
genetic cluster as the California locations. The assignment tests 
allowing computation of probabilities were run using Monte-Carlo 
resampling with N = 10,000 and a = 0.05 [20]. Results from the 
individual probability assignment tests were compiled and plotted 
using the software R v. 3.0.1. [12]. Additionally, individual 
assignment tests of 5 individuals from each of the reference 
populations and group self-assignment tests were performed to 
evaluate the accuracy of the assignment method. Group self- 
assignment tests were 100% accurate for all reference populations 
and genetic clusters. The individual approach considering the 
probabilities of computation assigned the highest probability to the 
correct population of origin for 64.5% of the individuals, while the 
correct genetic cluster was that of the highest probability for 
90.6% of the individuals. 

To test for possible sample size effects on STRUCTURE, we 
evaluated the possibility that uneven sample size might have an 
effect on the STRUCTURE analysis and our conclusions. We 
performed resampling of the data from all 20 North American 
samples, with and without replacement, with all samples equal to 
15 (without replacement) and 30 (with replacement) genotypes. 
When the genotypes were subsampled without replacement, 6 out 
of 10 runs resulted in K = 3 being the optimal number of clusters, 
while the remaining 4 had an optimal K = 2. When the genotypes 
were resampled with replacement, the optimal number of clusters 
was K = 3 for 7 out of 10 runs, K = 5 for 2 out of 1 0 runs, and 
K = 2 for 1 out of 1 0 runs. In all simulations, the California 
samples had greatest affinity to the Southeast US group samples at 
K = 3 (data not shown). 

Raw data (allele frequencies) can be found in Supplementary 
Materials, Table S4. 

Results 

Genetic diversity in Ae. aegypti from California is similar 
to that in natural populations 

Allele richness and heterozygosity values were estimated for 
both laboratory colonies and natural populations of Ae. aegypti, to 
address whether the California introduction was a consequence of 
a laboratory release (Table 1 and S3). The mean observed average 
heterozygosity (H„) for natural populations was significantly higher 
than that of laboratory strains (0.5153±0.076 and 0.3015±0.184 
respectively; P<0.001). Allelic richness across loci within natural 
strains was 3.6206±0.5681 and 2. 1633±0.6104 within lab strains, 
consistent with the genetic diversity pattern observed from the 
heterozygosity values. When the California populations were 
analyzed independendy, the H 0 values of all three California 
collections fell within the values of the natural populations (San 
Mateo: 0.4635; Madera: 0.5563; and Fresno: 0.4904; P = 0.7052; 
Table 1 and S3). 



PLOS Neglected Tropical Diseases | www.plosntds.org 



2 



July 2014 | Volume 8 | Issue 7 | e3029 



Aedes aegypti in California 



A) 



M 





gF$F $F £F & & & & & & ^ 



^> & ^> i> c> -i -b o" cr o ^ ^ ^ <? cr 



<o" <a <o <c 




B) 



S." O- O ^ <j. v., 



Q 



-5 




cT 



cT 



cT 



# ### # 



„cr 



/ 



j£ / ^ 



/ / / / 

/ / <f / / .# / / 



C? 



c) 




Figure 1 . Genetic structure within pantropical populations of Aedes aegypti. STRUCTURE bar plots indicating relatedness of Aedes aegypti 
populations based on 12 microsatellite loci. Each vertical bar represents an individual. The height of each bar represents the probability of 
assignment to each of K optimal clusters (different colors) determined using the Delta K method. (A) North America and Asian populations (K = 2), 
and (B) North American populations (K = 3). (C) Map indicating the North American geographic locations sampled in this study. (A) California, (O) 
other locations in North America. 
doi:1 0.1 371 /journal.pntd.0003029.g001 



California Ae. aegypti is closely related to populations in 
the southeastern USA 

To identify the potential source of the California introduction, 
we analyzed the pattern of population structure in a dataset that 
included the three California populations (Madera, Fresno, and 
San Mateo Counties) and a set of reference populations from the 
Americas and Asia (Table SI). Because California populations 
belong to the pantropical unit outside Africa known as Ae. aegypti 
aegypti [6, Fig. SI], in our comparisons we excluded Africa as a 
possible origin. Bayesian clustering analysis on these collections 



identified two main clusters (A K = 2) that split North America 
(USA and Mexico) from a group comprising the Caribbean, South 
America, and Asia (Figure 1A). Some admixture was observed 
between the groups, especially within the Florida populations. An 
independent Bayesian clustering analysis on the North American 
populations identified three clusters (A K = 3; Figure IB), corre- 
sponding with 1) Southern Mexico, 2) Arizona and Northern 
Mexico, and 3) southeastern USA, California, and the population 
in central Mexico (Amacuzac, Morelos). Therefore, California 
populations are genetically distinct from Ae. aegypti from Asia, 
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Table 1. Genetic diversity of Aedes aegypti populations. 





Locality 




H„ 


He 


AR 


San Mateo 




0.4635 


0.4949 


3.65 


Madera 




0.5563 


0.5435 


3.53 


Fresno 




0.4904 


0.5254 


3.67 


Pantropical* 


± SD 


0.5153±0.0760 


0.5 166 ±0.0764 


3.6206±0.5681 


Lab strains** 


± SD 


0.3015±0.1840 


0.3052±0.1963 


2.1633±0.6104 



H Q = observed heterozygosity; H e = expected heterozygosity; AR = Allelic richness estimated by rarefaction (N = 30 genes). 
*Pantropical = mean across populations from Asia and the Americas. 

**Lab strains = mean across Hamburg, Rockefeller, and Liverpool laboratory strains provided by David Severson (University of Notre Dame, Indiana). 
doi:1 0.1 371 /journal.pntd.0003029.t001 
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Figure 2. Individual and group mosquito genetic assignments. Percentage of individuals from Madera (A), Fresno (B), and San Mateo (C) 
counties assigned with the highest probability to each of the reference populations. (D) Scores calculated for each of the reference populations after group 
assignment of each of three California populations. Assignments were performed using Bayesian criteria for likelihood estimation with GENECLASS 2.0. 
doi:1 0.1 371 /journal.pntd.0003029.g002 
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Australia, Tahiti, Hawaii, and South America and thus these 
locales are highly unlikely to have been the source of the invasion. 
Interestingly, Ae. aegypli from California is most closely related to 
USA populations east of Arizona, rather than to Arizona 
populations, the state immediately adjacent to California 
(Figure IB). Geographic clustering can be observed in Figure 1C 
when we exclude the California populations. 

Likelihood estimation methods [18] used to more finely detail 
the origin of the California introduction, consistendy assigned 
California individuals to New Orleans with the highest probability 
(individual assignment tests; Fig. 2A-C); or California populations 
to Houston, with New Orleans the second most probable source 
(group assignment tests; Fig. 2D). 

In addition, we ran a Principal Components Analysis (PCA) on 
the allele frequencies of the 20 North American samples. The 
California samples are indistinguishable from samples from the 
Southeast US, but distinct from all other samples (Figure S2). 

Discussion 

Our analysis suggests that Ae. aegypli was introduced to 
California from the southeastern USA by multiple founder 
mosquitoes. Populations of recently introduced species are often 
relatively low in genetic variation due to the new populations 
having been started by a few founders. In comparing the level of 
genetic variation of California Ae. aegypti populations to other 
populations there is no indication of a decrease in genetic diversity 
(Table 1 and S3), including no decrease in the numbers of alleles at 
loci, the genetic parameter most sensitive to founder effects. The 
high level of genetic diversity observed in the California 
populations also argues against the release of a laboratory strain 
(Table 1 and S3). Additionally, California laboratories maintaining 
Ae. aegypti indicated that no recent strains of Ae. aegypti from the 
southeast USA are being currently reared in California, thus 
eliminating this alternative hypothesis. 

It is conceivable that the California populations are not recent 
introductions, which might also explain the relatively high level of 
genetic variation. We consider this less likely due to the fact that 
Ae. aegypti are morphologically quite distinctive and would have 
been detected given the extent of mosquito monitoring activities in 
California. The fact that sporadic Ae. aegypti in or near airports in 
California have been previously reported indicates that the quality 
and intensity of monitoring has been sufficient to detect this 
species. 

Our data indicates that the origin of the California invasion can 
be tracked to the region of New Orleans and Houston, cities 
~300 km apart, and likely did not originate from spillover from 
geographically adjacent regions that have long had Ae. aegypti 
(Arizona and northern Mexico). This suggests that commerce via 
air, railroad, or trucks from the New Orleans/Houston areas may 
have been responsible, though the mode of introduction remains 
unknown. 

The extensive mosquito surveillance program in California 
allowed for a rapid response to the initial detections, emphasizing 
the importance of enhanced surveillance and personnel training, 
especially in regions susceptible to invasion of exotic disease 
vectors. In the USA, dengue became a nationally notifiable 
infectious disease in 2010, and for 2013, California has reported 
1 24 imported cases [2 1] . In contrast to our previous study of an 
Ae. aegypti invasion into temperate Netherlands [22], it is 
conceivable that the climate of California may allow for these 
newly detected populations to overwinter and become established. 
The combination of imported dengue cases and the permanent 
presence of the primary vector of dengue in California could 



eventually lead to local transmission of the disease, thus becoming 
an important threat to public health. Identifying the origin of the 
California introduction should facilitate the implementation of 
strategies to reduce or eliminate the movement of new vectors into 
the area. 

The degree of accuracy in the identification of the source 
population of the California Ae. aegypti is proportional to the 
completeness of the reference population panel. Prior to this study, 
we had very sparse sampling from Arizona and northern Mexico; 
this delayed our analyses by several months while additional 
samples from this region were obtained and genotyped. Had our 
genetic panel been more complete, the study could have been 
carried out within a week. This highlights the importance of 
having in place extensive population genetic databases for invasive 
vectors prior to new invasions, which occur on a regular basis and 
are expected to increase as the climate patterns change. While 
spread of invasive vectors usually occurs inadvertendy, purposeful 
introductions by bioterror activities could be tracked in a similar 
manner. An analogy can be made with the notorious 200 1 anthrax 
event, whereby the origin of the strain was quickly traced thanks to 
an extensive database of DNA sequences of hundreds of anthrax 
genotypes [23]. More extensive reference databases of genetically 
fingerprinted human disease vector populations would greatly aid 
public health efforts to prevent the introduction and spread of 
vector-borne diseases. 

Supporting Information 

Figure SI California populations belong to the pantropical 
cluster of Aedes aegypti. STRUCTURE bar plots indicating 
relatedness of Aedes aegypti populations worldwide, based on 12 
microsatellite loci. Each vertical bar represents an individual. The 
height of each bar represents the probability of assignment to 
K = 2 optimal clusters (different colors) determined using the Delta 
K method. Pantropical cluster is shown in pink and the African 
cluster in blue. Populations from the African cluster were 
originally analyzed in Brown et al. 201 1 [6]. USA: United States 
of America, MEX: Mexico, COL: Colombia, BRA: Brazil, SA: 
Saudi Arabia, THA: Thailand, AU: Australia, CM: Cameroon, 
UG: Uganda, and SN: Senegal. 
(EPS) 

Figure S2 Principal Component Analysis in 20 North American 
samples of Aedes aegypti using 12 microsatellite loci. Note that 
California samples (*) are indistinguishable from samples from the 
Southeast US. 
(EPS) 

Table SI Collection information for Aedes aegypti populations 

analyzed for 12 microsatellite loci. 

(DOCX) 

Table S2 Number of alleles per microsatellite locus. 
(DOCX) 

Table S3 Genetic diversity of Aedes aegypti populations. 
(DOCX) 

Table S4 Raw allele frequencies at 12 microsatellite loci. 
(DOCX) 

Text SI Marker validation. 
(DOCX) 

Acknowledgments 

We thank collaborating mosquito and vector control agency staff who 
provided California samples: L. Irby, M. Garcia and A. Scalzo (Madera), 



PLOS Neglected Tropical Diseases | www.plosntds.org 



5 



July 2014 | Volume 8 | Issue 7 | e3029 



Aedes aegypti in California 



S. Mulligan and J. Holeman (Fresno), and R. Gay, N. Zahiri and B. Weber 
(San Mateo). Collections from outside of California not previously reported 
were kindly provided by B. A. Soria (Tijuana), M. Riehle (Hermosillo, 
Nogales), C. Gonzalez and the CERECOVE Panchimalco, SSM 
(Amacuzac), A. D. Lopez, A. D. Rodriguez, and the Vector Control 
Program of Chiapas (Mazatan and Tapachula), K. Smith (Maricopa, AZ), 

C. Ocampo (Cali), A. Malavasi (Jacobina), E. Khater (Saudi Arabia), A. 
Ponlawat (Bangkok), and D. Severson (lab strains). We thank B. Evans and 

D. Edwards for helpful discussions. 



Author Contributions 

Conceived and designed the experiments: AGS JRP. Performed the 
experiments: AGS JEB VK MHY. Analyzed the data: AGS JRP. 
Contributed reagents/materials/analysis tools: VK MHY. Wrote the 
paper: AGS JRP. 



References 

1. Aflez G, Rio M (2013). Dengue in the United States of Amcriea: A Worsening 
Scenario? BioMcd. Res. Int. 2013. ID678645. 1.3 p. dor.10.1 155/2013/678645 

2. Bhatt S, Gcthing PW, Brady OJ, Messina JP, Farlow AW ct al. (2013).The 
global distribution and burden of dengue. Nature 496: 504-507. doi: 10.1038/ 
nature 12060 

3. Lounibos LP (2002). Invasions by insect vectors of human disease. Ann. Rev. 
Entomol. 47: 233-266. 

4. Centers for Disease Control and Prevention Maps (2014). Available: www.cdc. 
gov/chikungunya/pdfs/CHIKV_VeetorControl. Accessed 8 January 2014. 

5. Jewell D, Grodhaus G (1984). Commerce and the Spread of Pests and Disease 
Vectors. Laird, M. Ed., New York: Prager. pp 103-107. 

6. Brown JE, McBridc CS, Johnson P, Ritchie S, Paupy C et al. (201 1). Worldwide 
patterns of genetic differentiation imply multiple 'domestications' of Aedes 
aegypti, a major vector of human diseases. Proc. Roy . Soc. Lond. B. 278: 2446- 
2454. doi: 10.1098/rspb.2010.2469 

7. Brown JE, Evans BR, Zheng W, Obas V, Barrera-Martinez L, et al. (2014). 
Human impacts have shaped historical and recent evolution in Aedes aegypti, 
the dengue and yellow Fever mosquito. Evolution. 68(2): 514-25. doi: 10.1111/ 
cvo.12281 

8. Raymond M, Rousset F (1995). GENEPOP (Version 1.2): Population Genetics 
Software for Exact Tests and Ecumenieism. J. Hered. 86: 248-249. 

9. Rousset F (2008). Genepop'007: a complete rc-implementation of the genepop 
software for Windows and Linux. Mol. Ecol. Resour. 8: 103-106. doi: 10.11 1 1/ 
j. 1 47 1-8286.2007.0 193 1 .x. 

10. Mcirmans PG, Van Ticnderen PH. (2004). GENOTYPE and GENODIVE: 
two programs for the analysis of genetic diversity of asexual organisms. Mol. 
Ecol. Notes 4: 792-794. DOI: 10. 1 1 1 1/j. 1471-8286.2004.00770.X 

1 1. Jombart T. (2008). Adcgenet: a R package for the multivariate analysis of genetic 
markers. Bioinformatics. 24: 1403-1405. doi: 10.1093/bioinformatics/btnl29. 



12. R Core Team. (2013). ISBN 3-900051-07-0. http://www.R-project.org/ 

13. Kalinowski ST. (2005). HP-RARE 1.0: a computer program for performing 
rarefaction on measures of allelic richness. Mol. Ecol. Notes. 5: 187-189. DOI: 
10. 1 1 1 1 /j. 147 1 -8286.2004.00845.x 

14. PritchardJK, Stephens M, Donnelly P. (2000). Inference of population structure 
using multilocus genotype data. Genetics 155: 945-959. 

15. Earl SA, vonHolt BM. (2012). STRUCTURE HARVESTER: a website and 
program for visualizing STRUCTURE output and implementing the Evanno 
method. Conserv. Genet. Resour. 4(2): 359-361. 

16. Evanno G, Regnaut S, Goudet J. (2005). Detecting the number of clusters of 
individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14: 
2611-2620. 

17. Rosenberg NA. (2004). DISTRUCT: a program for the graphical display of 
population structure Mol. Ecol. Notes. 4: 137-138. DOI: 10. 1 046/j. 147 1- 
8286.2003.00566.x 

18. Piry S, Alapctitc A, Cornuct JM, Pactkau D, Baudouin L, ct al. (2004). 
GENECLASS2: a software for genetic assignment and first-generation migrant 
detection. J. Hered. 95: 536-539. 

19. Rannala B, Mountain JL. (1997). Detecting immigration by using multilocus 
genotypes. Pror. Natl. Acad. Sci. U.S.A. 94: 9197-9201. 

20. Pactkau D, Slade R, Burden M, Estoup A. (2004). Genetic assignment methods 
for the direct, real-time estimation of migration rate: a simulation-based 
exploration of accuracy and power. Mol. Ecol. 13: 55-65. 

21. Centers for Disease Control and Prevention Maps (2013). Available: http:// 
diseascmaps.usgs.gov/dcp_us_human.html. Accessed 7 February 2014. 

22. Brown JE, Scholte E.-J., Dik M, Hartog WD, BeeuwkesJ, ct al (2011) Emerg. 
Infect. Dis. 17: 2335-2337. 

2.3. Hughes JM, Gcrberding JL. (2002). Anthrax bioterrorism: lessons learned and 
future directions. Emcrg. Infect. Dis. 8: 1013-2014. 



PLOS Neglected Tropical Diseases | www.plosntds.org 



6 



July 2014 | Volume 8 | Issue 7 | e3029 



